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Abstract 

Background: Measurements of neutron radii provide important constraints on the isovector 
sector of nuclear density functionals and offer vital guidance in areas as diverse as atomic parity 
violation, heavy-ion collisions, and neutron-star structure. Purpose: To assess the impact of spin- 
orbit currents on the electromagnetic- and weak-charge radii of a variety of nuclei. Special emphasis 
is placed on the experimentally accessible electroweak skin, defined as the difference between weak- 
charge and electromagnetic-charge radii. Methods: Two accurately calibrated relativistic mean 
field models are used to compute proton, neutron, charge, and weak-charge radii of a variety of 
nuclei. Results: We find that spin-orbit contributions to the electroweak skin of light neutron-rich 
nuclei, such as 22 and 48 Ca, are significant and result in a substantial increase in the size of the 
electroweak skin relative to the neutron skin. Conclusions: Given that spin-orbit contributions 
to both the charge and weak-charge radii of nuclei are often as large as present or anticipated 
experimental error bars, future calculations must incorporate spin-orbit currents in the calculation 
of electroweak form factors. 
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I. INTRODUCTION 



Recently the Lead Radius Experiment ( "PREX" ) at the Jefferson Laboratory used par- 
ity violating electron scattering to probe the weak-charge density of 208 Pb [TJ |2]. In Born 
approximation, the parity violating asymmetry — the fractional difference in cross sections 
for positive and negative helicity electrons — is directly proportional to the weak form factor, 
which is the Fourier transform of the weak-charge density. Although for a heavy nucleus one 
must include the effects of Coulomb distortions, these have been accurately calculated [3j H] . 
Thus, parity-violating electron scattering can be used with as much confidence to measure 
the weak-charge form factor of the nucleus as conventional electron scattering has been 
used throughout the years to accurately map the electromagnetic charge distribution. Many 
details of a practical parity- violating experiment, along with a number of theoretical correc- 
tions, have been discussed in Ref. [5|. PREX demonstrated excellent control of systematic 
errors and showed that accurate and model independent measurements of weak-charge den- 
sities are now feasible. 

In this paper we focus on the electroweak skin of a nucleus which we define as 

-Rwskin = Rwk Rch j ( 1 ) 

where i? w k and R c ^ are the root-mean-square radii of the weak-charge and electromagnetic- 
charge densities, respectively, Note that unlike the point-neutron and point-proton densities, 
the weak-charge density — the source for the Z° weak boson — and the electromagnetic-charge 
density are physical observables. For simplicity, we refer to -R W skin as the weak skin. Given 
that the charge radius of 208 Pb is very accurately known [6], the determination of the weak 
form factor of 208 Pb by the PREX collaboration resulted in the following value for its weak 
skin [2] : 

#wskin( 208 Pb) = 0.32 ± 0.18(exp) ± 0.03(mod) fm . (2) 

Note that the experimental error ( "exp" ) includes both statistical and systematic errors while 
the small model error ( "mod" ) describes the sensitivity in the extraction of i? w k due to un- 
certainties in the surface thickness [2] . Although the statistical accuracy of the measurement 
was compromised by unforeseen technical difficulties, a follow-up measurement ( "PREX-II" ) 
designed to achieve the original ±0.05 fm goal has been proposed and accepted [7J. Moreover, 
a fresh new proposal has been made to use parity-violating electron scattering to measure 
the weak form factor of 48 Ca at a momentum transfer of q = 0.778 fm -1 . The Calcium 
Radius Experiment ("C-REX") is designed to constrain the weak-charge radius of Ca to 
an accuracy of ±0.03 fm [8J. A measurement in a smaller neutron-rich nucleus is desirable 
because the form factor can be measured at a larger momentum transfer where the parity- 
violating asymmetry is larger; for comparison, the weak form factor in 208 Pb was measured 
at q = 0.475 fm -1 . In addition, a precise measurement of the weak radius of 48 Ca may 
have a significant impact on nuclear structure as it provides information that is indepen- 
dent and complementary to that found in 208 Pb [U] . Furthermore, it is reasonable to expect 
first-principle calculations of the structure of 48 Ca where the role of 3-neutron forces may 
be particularly interesting and important. Finally, in conjunction with PREX and PREX- 
II, C-REX will provide vital guidance in areas as diverse as atomic parity violation [JUl ITT] , 
heavy-ion collisions [T2HT6] . and neutron-star structure [T7J422]- Given the expected accuracy 
of these pioneering experiments, it is critical to assess the role of u sub-leading" contributions 
to the weak skin of these nuclei. Thus, it is the main goal of the present contribution to 
quantify the impact of spin-orbit currents on the electroweak skin of a variety of nuclei. 
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We note that meson-exchange currents — which are not considered in this contribution — can 
change the distribution of both electromagnetic and weak charge. For example pn^ and 
i>jo^( meson-exchange currents (MEC) are known to modify significantly the electromagnetic 
form factors of the deuteron [23] • However, heavy mesons are unlikely to transport charge 
over large distances, so we do not expect significant MEC corrections to either R ch or i? wk 
for a heavy nucleus. Yet MEC corrections could be more important for the case of the weak 
skin (i? w k — -Rch) and this deserves further consideration. 

The result depicted in Eq. ^ represents a true experimental milestone. It provides direct 
experimental evidence that the weak-charge density in 208 Pb extends further out than the 
corresponding electromagnetic-charge density. That is, there is an enhancement of weak 
charges (which are dominated by neutrons) relative to electromagnetic charges (which are 
dominated by protons) near the nuclear surface. A quantity that is closely related to the 
weak skin of a heavy nucleus is the neutron skin. In analogy to the weak skin, the neutron 
skin is defined as the difference between the point neutron R n and point proton R p root- 
mean-square radii: 

-Rnskin = Rn Rp ■ (3) 

The weak and neutron skins are closely related to each other because the weak charge of 
a neutron (« —1) is much larger than the weak charge of the proton (« 0.07) largely in 
the same way as the electromagnetic charge of the proton is much larger than that of the 
neutron. Indeed, to a good approximation the weak form factor of a single nucleon equals 
the negative of the electromagnetic form factor of its isospin partner (see Appendix). The 
structure of a nucleus is often modeled in terms of point nucleons interacting via strong 
nuclear potentials or meson exchanges. Therefore, nuclear-structure calculations provide 
predictions for the neutron skin whereas it is the weak skin that is experimentally accesible. 
In the present paper we calculate the difference between iCskin and -Rnskin that originates 
from the internal structure of the nucleon. The internal structure of the nucleon is contained 
in a variety of well-measured single nucleon electromagnetic and weak form factors. The 
electromagnetic size of the nucleon leads to a charge radius of the nucleus i? c h that is slightly 
larger than the point proton radius R p due (mostly) to the charge radius of a single proton r p 
(-R^ h ~-Rp + rp). Similarly, the weak radius _R w k is slightly larger than R n (mostly) because of 
the weak-charge radius of a single neutron. Note that we have reserved the use of upper-case 
R to nuclear radii and of lower-case r to the radii of single nucleons. 

The manuscript has been organized as follows. In Sec. [IT] we outline the formalism used 
to calculate weak and electromagnetic charge form factors and densities. These form factors 
are computed from one-body currents that include weak and electromagnetic single-nucleon 
form factors. We discuss in detail the contribution to both of the form factors from the 
spin-orbit currents. For a recent reference on the impact of spin-orbit currents on the charge 
density of light nuclei see Ref . [21] . However, to our knowledge the role of spin-orbit currents 



on the weak-charge density has never been studied. In Sec. Ill we present results for proton, 
neutron, charge, and weak radii for a variety of nuclei. Spin-orbit corrections are particularly 
large for light nuclei and we observe a significant increase in the weak skin of these nuclei 
relative to their neutron skin. We offer our conclusion in Sec. 
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II. FORMALISM 



We start the formalism by writing the most general form of the single-nucleon matrix 
elements of the electroweak current consistent with Lorentz covariance and parity invariance. 
That is I 



(N(p' y s')\J£ M \N(p,s)) = U(p' y s') 
(N(p>,s')\j£ c \N(p,s)) = U(p',s') 



F^ + iF 2 a- 



2M 
2M. 



U(p,s) 
U(p,s) 



(4a) 
(4b) 



where U(p,s) are on-shell nucleon spinors, q = p'—p is the four momentum transfer to 
the nucleon, and ^1,2(^1.2) are electromagnetic (weak-neutral vector) Dirac and Pauli form 
factors, respectively Note that the axial-vector component of the weak neutral current 
makes no contribution to the elastic form factor of a nuclear ground state of definite parity. 
As can be seen from the above equations, the formal structure of the vector current is 
the same regardless of whether one considers electromagnetic or weak matrix elements or 
whether one considers protons or neutrons. Thus, for illustration purposes we focus the 
derivation on the proton contribution to the charge form factor. 



A. Charge Form Factor 



Throughout this contribution we will assume the validity of the impulse approximation, 
namely, that the single-nucleon form factors defined in Eq. ^ may be used without mod- 
ification in the nuclear medium. Moreover, we rely on an alternative description of the 
single-nucleon electromagnetic current based on the electric and magnetic Sachs form fac- 
tors rather than on F\ 2. The two sets of form factors are related as follows: 



G E (Q 2 ) = F^Q 2 ) - tF 2 (Q 2 ) 
G m (Q 2 ) = F^Q 2 ) + F 2 (Q 2 ) 



(5a) 
(5b) 



where Q 2 = —q^ >0 and r = Q 2 /AM 2 . In turn, the electromagnetic current operator may 
be written as 



Gm — Gi 

1+T 



T7 M + ia 



fJ,U 



2M 



(6) 



Given that we are interested in computing the elastic form factor of spherical nuclei, only 
the density (i.e., the zeroth component of the current) needs to be consider. Such density is 
given by 



90 

J EM 



g e7 + 



1 + T 



. ' ' 2M 



(7) 



Thus, in the impulse approximation, the charge form factor of the nucleus may written as 
follows: 



ZF ch (q) = G E (q 2 )F y (q) + 



G M (q 2 )-G E (q 2 

1+T 



rF Y (q) + ^F T (q) 



(8) 



where Z is the number of protons and the charge form factor of the nucleus has been 
normalized to F c h(q = 0) = 1. Note that all nuclear structure information is contained in the 
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elastic (point-nucleon) vector and tensor form factors defined as 

F y (q) = J ^(r)e iq - r 7 V(r)rf 3 ^ , (9a) 

F T (q) = J ^(r)e iq ' r 7°7 • q^(r)rf 3 r . (9b) 

We now proceed to evaluate these two nuclear form factors using a relativistic mean-field 
approximation. 

B. Relativistic Mean-Field Approximation 

In a relativistic mean-field approximation, the eigenstates of the Dirac equation corre- 
sponding to a spherically symmetric ground state may be classified according to a generalized 
angular momentum k [26]. That is, the single-particle solutions of the Dirac equation may 
be written as 

K-.w=zC-H:Mn . do) 



r V ifnK(r)y- Km (r) 

where n and m are the principal and magnetic quantum numbers, respectively, and the 
spin-spherical harmonics are defined as follows: 

_ . /~\ / A I J • \ ■ II 1 7 & J if/x^O^ f-tt\ 

y Km {r) = (Y\l\jm) ; j = \k\--\ I = < _^ _ ^ if ^ < Q ^ ' 



Note that for the phase convention adopted in Eq. ( 10 ) (namely, the relative factor of %) 
both g{r) and f(r) are real functions in the case of real mean-field potentials. Further, the 
following spinor normalization has been used: 

Ui Km {v)U nKm {v)d?r = [g 2 nK (r) + fl(r)\dr = 1 . (12) 

Thus, in a mean-field approximation the vector and tensor form factors may be expressed 
as a sum over occupied single-particle states. That is, 

F v(q) = / W n K m(r)j (gr)7°W nKm (r)d 3 r , (13a) 

num. 

Ft(q) = J2 I UnKmirKiqr) *7°7 • rW nKm (r)c/ 3 r . (13b) 

nnm 

Here j l (qr) are spherical Bessel functions and the sum over single-particle quantum numbers 
{num} is restricted to the occupied orbitals. Note that by virtue of the spherical symmetry of 
the nuclear form factors, all explicit dependence on q has been eliminated through an "angle- 
average" (i.e., by integrating over ciq/47r). This fact leads to a considerable simplification 
and enables one to readily perform the integrals over the remaining solid angle. One finally 
obtains: 

POO 

F y (q) = X)(2J + 1) / dL(r) + fl(r) j (qr)dr , (14a) 

/■oo 

F T (q) = y>(2j+l) / g nK (r)f nK (r)jAqr)dr. (14b) 

Jo 
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Note that for a spherically symmetric ground state there is a third independent "scalar" 
form factor that is identical to the vector one except for a relative minus sign between upper 
and lower components. Given the form of the electromagnetic current adopted in Eq. Q, 
and ignoring possible off-shell ambiguities, the scalar form factor plays no role in the present 
discussion. 

In summary, nuclear charge and weak form factors — both normalized to one at q — — are 
now obtained from Eq. pi) by properly adding proton and neutron contributions. That is, 



ZF ch (q) 



i=p,n 




G^q 2 )F^q) + 
G^q 2 )F^q) + 



i=p,n 



'Gk(q 2 )-Gi{q 



GW)-G£(g 
1 + r 




TF4M + ^F4(q) 



TF4(q) + ^Fi(q) 




where Q wk is the weak- vector charge of the nucleus (see Appendix). 



C. Electromagnetic and Weak Charge Radii 



The experimentally measurable electromagnetic and weak i? w k charge radii are ob- 
tained from the slope of the respective form factors at the origin. That is, 



R 



ch 



-6 



dF, 



cli 



dq 2 



and R 



q=0 



wk 



-6- 



dF 



wk 



dq 2 



q=0 



(16) 



Although in the next section we provide results in the form of tables and figures, it is 
illuminating to di scuss the particle content of both i? c h and -R w k- We start with the former. 
According to Eq. (15a), the electromagnetic charge form factor to first order in q 2 may be 
written as follows: 



6 

+ N [ 1 - q ^Rl 



ZF ch (q) = Z ( 1 - ^Rl ) ( 1 - £rj ) + rb^(l + ft) + • • • 



6 
n 2 



K + 



AM 2 

±_ 
AM 2 



n n N{l + P) + 



(17) 



where r 2 n and fcp } Ti 3X6 Sill gle-nucleon electromagnetic mean-square radii and anomalous 
magnetic moments, respectively. Note that we have defined 



J T 



4M 
AM 



E(^' +1 ) / 

Jo 



rgL(r)fL(r)dr . 



rgL(r)f: K (r)dr. 



(18a) 
(18b) 



of the nucleus. That is, 



Using Eq. (|17|) one can readily obtained an expression for the electromagnetic charge radius 

(19) 
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where "spin-orbit" contributions to the charge radius have been defined as follows: 



(^)so = -^(l + / T p ), (20a) 



(^)so = -|^(l + / T n ). (20b) 

In the context of a relativistic mean-field approximation — where both upper and lower 
components are self-consistently generated from the Hartree equations — there are no further 
simplifications. However, one can shed light into the nature of the spin-orbit contribution 
by generating the lower component from the upper component by assuming a free-space 
relation. That is, 

f ^ = m(Tr + v)^ r) - (21) 

In this free-space limit an enormous simplification ensues, as the tensor integrals given in 



Eq. ( 18 ) can be evaluated in closed form. In this limit the spin-orbit contributions to the 



charge radius become 

«>» = -^(! + £) £ + 1) , (22a) 

riK 

(rlU = + /") £sgnW + 1) . (22b) 

nn 

The above results indicate that: (a) there is no contribution from s-states and (b) there is 
an exact cancellation between spin-orbit partners, which have the same I but opposite signs 
for k. Thus, in the particular case of 48 Ca the only spin-orbit contribution comes from the 

"unpaired" H7/2 neutron orbital. For this case we obtain the following figure of merit: 

(Ac -f(Oso = |^- -0.101 fm 2 . (23) 

Thus, the spin-orbit contribution amounts to approximately — 0.015 fm of the total charge 
radius of 48 Ca, which is significantly larger than the quoted experimental error [5J. Note 
that in a self-consistent RMF approximation the cancellation between spin-orbit partners, 
although still large, will be incomplete. However, there is an additional cancellation between 
neutron and proton orbitals with the same quantum numbers due to their almost equal but 
opposite anomalous magnetic moments. Ultimately, the spin-orbit contribution is dominated 
by the unpaired orbitals. 

To obtain the weak-charge radius R wk one proceeds in an analogous manner but now 



starting from Eq. (15b). Expanding the weak charge form factor to first order in q one 
obtains 

Q^F wk (q) = Z (l- q jR^j gl (l - p^j + ^*pZ(l + f P T ) + ... 

+ * (l - jRl) £ (l " jrl) + ^p^N(l + £) + ... (24) 

where g% = 0.0721 and = —0.9878 are (radiatively-corrected) single-nucleon weak-vector 
charges 0, [25] . Note that the weak-charge radii r 2 p n and anomalous weak-magnetic moments 
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^p,n 



may be expressed exclusively in terms of the corresponding electromagnetic and strange- 
quark quantities, as indicated in the Appendix. In this way, the weak-charge radius of the 

. T 

N 



nucleus can now be readily extracted from Eq. (24). That is 

Z 



Q 



wk 



(K + fj) + <fj>. 



+ 



Q 



wk 



9 n v (Rl + fl) + iO* 



(25) 



where spin-orbit contributions to the weak-charge radius are defined in analogy to Eq. (22): 



OK r , 



2M 2 
~2NP 



(26a) 
(26b) 



Given that Q w ^ = Zg^+Ng^ (with « and ~ — 1) the weak-charge radius of the nucleus 
i? w k is dominated by the contribution from the neutron radius R n . 

We close this section by providing a figure of merit for the weak-charge radius of 48 Ca 
by neglecting the strange-quark contribution and by assuming the free-space relation in the 



evaluation of the tensor integrals [Eq. (18)]. We obtain: 



P 2 



Q 



wk 



glR A p + g y p + g. 



,n„2 



+ 



N 



Q 



wk 
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g»R 2 n + 9 y p + girl + i^{g> P + 9> n ) 



(27) 



Note that in the above expression single-nucleon mean-square radii (r^ n ) as well as anoma- 
lous magnetic moments (K p , n ) are purely electromagnetic. Assuming (point) proton and 
neutron radii for 48 Ca as predicted by the FSUGold model [27j (or "FSU" for short) we 
obtain 

i? wk = 3.679(3.669) fm , (28) 

where the quantity in parenthesis represents the FSU prediction without including the spin- 
orbit contribution that amounts to: 



N 



r 2 ) 



Q 



wk 



N 6(g»K p 



Q 



wk 



7M 2 



+0.077 fm 2 



(29) 



Given that the weak spin-orbit contribution has the opposite sign as the corresponding 
electromagnetic one [see Eq. (23)], the spin-orbit contribution to the weak skin of 48 Ca is 
significant; of about 0.03 fm (see Table[ll|) . This is within the projected accuracy of the 
proposed C-REX experiment [S]. 



D. Electroweak Nucleon Form Factors 

Since our main objective is the calculation of electromagnetic- and weak-charge radii, we 
adopt a simple dipole parametrization of the single nucleon form factors that is accurate 
at moderate values of the momentum transfer. That is, we define electromagnetic single 
nucleon form factors as follows 128 1: 



G|(Q 2 ) = = = G D (Q 2 ) , (30) 

ftp ftn 
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where the dipole form factor is given by 

Gb(Q 2 ) = (l + ) . (31) 



Here = 0.769 fm is the mean-square proton radius, fi p = 2.793 the proton magnetic 
moment, and fi n — —1.913 the neutron magnetic moment. For the electromagnetic neutron 
form factor — which vanishes at Q 2 = — we rely on the following Galster parametrization 



G ^ 2 > - - {Swm) Gd(q2) ■ (32) 

where r 2 = — 0.116 fm 2 is the electromagnetic mean-square radius of the neutron. Note 
that all form factors are expressed exclusively in terms of experimentally determined single 
nucleon mean-square radii and magnetic moments [29J . Finally, given that the strange-quark 
form factors of the nucleon are small at small momentum transfers, we will ignore them in 
this contribution. In this case, the weak-charge form factors of the nucleon can be expressed 
exclusively in terms of the corresponding electromagnetic form factors given here and the 



weak charges of the nucleon (see Appendix). For example, using Eqs.(A.8) and (A. 9) the 
weak mean-square radii and magnetic moments are given by 

?l = 2.358 fm 2 , fl p = +2.091 ; (33a) 

f\ = 0.777 fm 2 , fl n = -2.897 . (33b) 



III. RESULTS 

In this section we present results — primarily proton, neutron, charge, and weak-charge 
radii — for a variety of nuclei as predicted by two accurately calibrated relativistic mean field 
models: FSU [27] and NL3 151] . We start by displaying in Table|l]the contribution from 
the individual single-particle orbitals in 48 Ca to the charge and weak-charge spin-orbit radius. 
For illustration purposes, the predictions have been made using only the FSU parametriza- 
tion. The quantities in parenthesis represent the results obtained by assuming a free-space 



relation between upper and lower components, as indicated in Eq. (21). Qualitatively, all 
major trends in the results may be understood using this simplified case. In particular, the 
results displayed in parenthesis in Tablejl} (i) are independent of the dynamics, (ii) scale 



as 1(1+1), and (iii) display an exact cancellation among spin-orbit partners [see Eqs. (22)]. 
Thus, in this limit the sole contribution to the spin-orbit radius comes from unpaired spin- 
orbit partners (an unpaired If 7/2 neutron orbital in the case of 48 Ca). Given the natural 
ordering of spin-orbit partners in nuclei, namely, the k < orbital more deeply bound than 
the k > orbital, the spin-orbit contribution to the charge radius is always negative in the 
case of unpaired neutrons. Note that the situation is reversed in the case of protons due to 
an anomalous magnetic moment of opposite sign to that of the neutron. Moreover, in the 
case of the weak-charge radius, the spin-orbit contribution is always of opposite sign — for 
both neutrons and protons — than in the electromagnetic case. Indeed, for a given nucleon 
orbital the ratio of spin-orbit contributions is given by the following simple expression: 



(^p,n) so ( K pn \ ( Z \ ( K np \ ( Z 



< . (34) 
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Hence, although in general small, charge and weak-charge spin-orbit radii contribute with 
the same sign to the weak skin of a nucleus. For example, as indicated in Tableful the 
spin-orbit contribution amounts to about 0.04 fm in the case of the weak skin of 22 and to 
about 0.03 fm for 48 Ca. 



nlj (k) 


/ 2\ 
\ r p/so 


/ 2 \ 
Vn/so 




1si/ 2 (-1) 
lPs/2 (-2) 

lPl/2 ( + 1) 

lds/2 (-3) 
ld 3 /2 (+2) 
2s 1/2 (-l) 
%2 (-4) 


+0.570(+0.000) 
+3.107(+1.583) 
-1.997(-1.583) 
+7.138(+4.750) 
-6.205(-4.750) 
+0.178(+0.000) 
+0.000(+0.000) 


-0.648(+0.000) 
-3.409(-1.690) 
+2.113(+1.690) 
-7.710(-5.069) 
+6.613(+5.069) 
-0.177(+0.000) 
-13.042(-10.138) 


-0.078(+0.000) 
-0.302(-0.107) 
+0.116(+0.107) 
-0.572(-0.319) 
+0.408(+0.319) 
+0.001(+0.000) 
-13.042(-10.138) 


Total 


+2.791(+0.000) 


-16.260(-10.138) 


-13.469(-10.138) 


nlj («) 


(^p)so 




(r 2 )so 


lsi/ 2 (-l) 

1P3/2 (-2) 
lPl/2 ( + 1) 

ld 5/2 (-3) 
ld 3 /2 (+2) 
2s 1/2 (-l) 
%2 (-4) 


-0.489(+0.000) 
-2.669(-1.360) 
+1.716(+1.360) 
-6.132(-4.081) 
+5.331(+4.081) 
-0.153(+0.000) 
+0.000(+0.000) 


+0.494(+0.000) 
+2.595(+1.286) 
-1.608(-1.286) 
+5.869(+3.859) 
-5.034(-3.859) 
+0.135(+0.000) 
+9.928(+7.717) 


+0.005(+0.000) 
-0.074(-0.074) 
+0.108(+0.074) 
-0.263(-0.222) 
+0.297(+0.222) 
-0.018(+0.000) 
+9.928(+7.717) 


Total 


-2.396(+0.000) 


+12.379(+7.717) 


+9.983(+7.717) 



TABLE I: Contributions from the individual single-particle orbitals to the spin-orbit component of 
the mean square charge radius (upper table) and weak-charge radius (lower table) of Ca. Mean 
square radii are expressed in units of 1CU 2 fm 2 and were generated using the FSU interaction. 
Quantities displayed in parenthesis have been computed using the free-space relation given in 
Eq.pil. 



Also shown in Table[T] are predictions from the self-consistent RMF approach, where 
now both upper and lower components are dynamically generated. The qualitative trends 
discussed previously are clearly preserved, although now the cancellation among spin-orbit 
partners is incomplete. For example, the contributions from the lp orbitals differ from 
each other by about 55-60%. However, although the spin-orbit cancellation is incomplete, 
there is an additional cancellation stemming from the nucleon anomalous magnetic moments 
(which are nearly equal in magnitude but opposite in sign). Ultimately, the spin-orbit radius 
continues to be dominated by the unpaired lf 7 / 2 neutron orbital. Note, however, that the 
exact RMF prediction exceeds by about 30% the analytic result obtained from assuming the 
free-space relation. This enhancement is due to the reduction of the effective nucleon mass 
in the nuclear medium. 

In Tableful predictions are displayed for the root-mean-square radii for a variety of magic 
(or semi-magic) nuclei as predicted by the NL3 and FSU models. Also shown (last three 
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ATT O 


4.194 


A or*o 

4.308 


4.254(4.268) 


/I AC\ A I A OO0\ 

4.404(4.393) 


f\ A A A 

0.114 


0.149(0.126) 




T7CTT 

b SU 


A 1 O 1 

4.181 


4.269 


4.242(4.255) 


4.364(4.353) 


n no o 

0.088 


n 1 oo/n nno\ 

0.123(0.098) 


118c 
°Sn 


ATT O 


4.561 


A 'i n c\ 

4.760 


4.636(4.628) 


4.835(4.843) 


0.199 


r\ -\ C\C\ i C\ 01 CT \ 

0.199(0.215) 






4.559 


A '1C\'1 

4.707 


4.634(4.625) 


/I '70n/ /l ^700\ 

4.780(4.788) 


r\ A A o 

0.148 


n i A'i ( r\ i ro\ 

0.147(0.162) 


Sn 


ATT O 


4.o43 


4.989 


4.700(4.705) 


5.077(5.074) 


0.346 


0.377(0.369) 




■POTT 

FSU 


A C f A 

4.654 


4.925 


4.710(4.716) 


5.011(5.008) 


0.271 


0.301(0.292) 


208p^ 


ATT O 


5.46L) 


5.740 


5.510(5.514) 


5.815(5.814) 


n oon 

0.280 


n onr /n onn\ 

0.305(0.300) 




■priTT 

FSU 


5.469 


5.676 


5.519(5.523) 


5.749(5.747) 


0.207 


0.230(0.224) 


138 13 
Ba 


ATT O 


4.77b 


5.012 


4.827(4.838) 


cr ac\c\I^ nr»o\ 

5.100(5.093) 


f\ ,) ■ ) "-7 

0.237 


0.273(0.255) 




b SU 


4.775 


4.957 


4.826(4.837) 


5.043(5.035) 


n 1 o o 

0.182 


0.217(0.198) 


158 Dy 


NL3 


5.039 


5.235 


5.099(5.098) 


5.309(5.311) 


0.196 


0.210(0.212) 




FSU 


5.027 


5.172 


5.087(5.086) 


5.245(5.246) 


0.146 


0.158(0.160) 


176 Yb 


NL3 


5.215 


5.497 


5.273(5.272) 


5.573(5.574) 


0.282 


0.300(0.302) 




FSU 


5.208 


5.424 


5.266(5.265) 


5.498(5.498) 


0.215 


0.232(0.233) 



TABLE II: Proton, neutron, charge, and weak-charge radii (in fm) of a variety of nuclei as predicted 
by the NL3 and FSU relativistic mean field models. The last two columns display the neutron skin 
(R n -R p ) and weak skin (Rwk-Rch), respectively. Quantities displayed in parenthesis have been 
computed without the spin-orbit correction. 

entries) are nuclei of relevance to the atomic-parity violating program [II]. These nuclei 
are relevant as they are members of long chains of naturally occurring isotopes that help 
eliminate uncertainties in the atomic theory by forming suitable ratios of parity-violating 
observables. By doing so, uncertainties in neutron radii become the limiting factor in the 
search of physics beyond the standard model. Note that the NL3 and FSU models — although 
both accurately calibrated — predict large differences in neutron radii. However, we caution 
that in the case of 158 Dy and 176 Yb, nuclear deformation and pairing correlations — which 
have not been included — may play an important role. 

The first three columns of numbers in Table|IT| provide predictions for the point-proton, 
point-neutron, and charge radii of several nuclei. To a large extent these results are expected 
and in some particular cases (such as 208 Pb) have been extensively discussed. For example, 
given that the charge form factor of many of these nuclei has been accurately measured 
via (parity conserving) elastic electron scattering, this information — mostly in the form of 
charge radii — has been incorporated into the calibration of the RMF models. Thus, the 
predictions of both models for charge radii are in fairly good agreement (e.g., they differ by 
only 0.009 fm in the case of 208 Pb). In contrast, the lack of reliable neutron form factors 
leaves the isovector sector of the relativistic functionals largely unconstrained, thereby gen- 
erating large differences in the predictions of neutron radii (e.g., 0.064 fm, or more than 1%, 
in the case of 208 Pb). Fortunately, the prospects of constraining the isovector sector through 
the measurement of neutron-radii are very good. Indeed, the vigorous and highly success- 
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ful parity-violating program at JLAB has recently provided the first model-independence 
evidence in favor of a weak skin in 208 Pb. Moreover, a follow-up measurement has been 
approved to achieve the original ±0.05 fm goal (PREX-II) and a fresh new one (C-REX) has 
been proposed to constrain the neutron radius of 48 Ca to ±0.03 fm. It is therefore critical 
and timely to assess the impact of the spin-orbit corrections on the weak radius of these 
nuclei — which until now has never been considered. 
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FIG. 1: (Color online). Charge and weak-charge densities for 22 O and 48 Ca as predicted by the 
FSU interaction. The dashed lines represent the corresponding densities in the absence of spin-orbit 
corrections. 



Given the enormous cancellation among spin-orbit partners and proton-neutron orbitals 
with the same quantum numbers, the spin-orbit contribution to the charge and weak-charge 
radii lacks the coherence displayed by the dominant vector form factor [see Eqs. (14)]. Thus, 
one expects that the spin-orbit correction will be largest for light nuclei with unpaired 
spin-orbit partners. Indeed, Table|TT] displays relatively large spin-orbit contributions — of 
~ 0.045 fm and 0.03 fm — to the weak skin of 22 and 48 Ca, respectively. These spin-orbit 
corrections are commensurate with the projected experimental uncertainty so they must be 
included in the prediction of both charge and weak-charge radii. Charge and weak-charge 
densities for 22 and 48 Ca — with and without spin-orbit corrections — as predicted by the 
FSU interaction are also shown in Fig.[TJ Note that the effect from neglecting the spin- 
orbit contribution is clearly discernible in the figure, as is the fact that the modification 
to the charge and weak-charge radii goes in the opposite direction — thereby enhancing its 
contribution to the weak skin. As alluded earlier, the impact of the spin-orbit contribution 
diminishes with increasing baryon number and amounts to only « 0.002 fm in the case 
of the weak radius of 208 Pb — significantly below the anticipated ±0.05 fm of the PREX-II 
measurement. Yet, the ~ 0.004 fm contribution to the charge radius of 208 Pb is significantly 
larger than the minute 0.0009 fm error quoted in Ref. [6]. In the particular case of atomic 
nuclei of relevance to the atomic parity violating program, the spin-orbit contribution to 
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the weak-charge radius of 138 Ba, 158 Dy, and 176 Yb amounts to ~ 0.008 fm, ^— 0.002 fm, and 
~ —0.001 fm, respectively. Note that whereas the NL3 and FSU predictions for the weak skin 
of the various nuclei differ significantly, the spin-orbit contribution appears to have very little 
model dependence. In general, we find that the intrinsic structure of the nucleon leads to a 
weak skin that is larger than the corresponding neutron skin. Moreover, for light neutron-rich 
nuclei — such as 22 and 48 Ca — the spin-orbit contribution generates a further enhancement 
of the weak skin. In particular, whereas the FSU interaction predicts a neutron skin in 22 of 
-Rnskm = 0.417 fm, the prediction for the experimentally measurable weak skin is significantly 
larger, namely, R ws km = 0.487 fm. Although not as large, the effect is still significant in 
48 Ca [R nski n = 0.197 fm and R wskin = 0.247 fm) and even in 208 Pb {R nskin = 0.207 fm and 
-Rwskm = 0.230 fm). We also display some of these results in graphical form in Fig.[2] where 
predictions for the weak skin — with and without spin-orbit contributions — are shown as 
a function of the corresponding neutron skin. Again, the neutron skin is interesting as 
it represents a pristine nuclear-structure observables that is instrumental in constraining 
the isovector sector of the nuclear density functional. The weak skin, on the other hand, 
although sensitive to the internal structure of the nucleon, is both experimentally accessible 
as well as strongly correlated to the neutron skin. Note that in Fig. [2] we have added the 
line .Rwskin = -Rnskin to indicate that in all cases (except for 118 Sn) the weak skin is larger 
than the neutron skin — and significantly larger for the lighter neutron-rich nuclei. Finally, 
the fact that the points are more "compressed" in the case of FSU than NL3 is a reflection 
of the softer symmetry energy of the former relative to the latter. 
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FIG. 2: (Color online). Electroweak skin (i? w k--Rch) with and without spin-orbit corrections as 
a function of neutron skin (R n -R p ) for the various neutron-rich nuclei considered in this work. 
Predictions are made using both the NL3 and FSU interactions. 
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IV. CONCLUSIONS 



A relativistic mean-field approximation has been used to compute proton, neutron, 
charge, and weak-charge densities and form factors of a variety of neutron-rich nuclei. Spe- 
cial emphasis has been placed on the impact of spin-orbit currents on the electroweak skin 
of these nuclei. Although closely related to the neutron skin, the weak skin — defined as the 
difference between the experimentally accesible weak- and electromagnetic-charge radii — is 
sensitive to the internal structure of the nucleon. The weak-charge radius of a nucleus is 
closely related to the (point) neutron radius because the weak charge of the neutron is much 
larger than that of the proton. This is analogous to the reason why the charge and pro- 
ton radii are closely related. In the absence of spin-orbit corrections, and regardless of the 
model, we found a weak skin larger than the neutron skin for all nuclei investigated here. 

Once spin-orbit currents were incorporated, a significant increase in the weak skin- 
especially in the case of the two lighter nuclei 22 and 48 Ca — was observed. In particular, 
we found quite generally that the spin-orbit contribution to the weak and charge radii 
enter with opposite sign so their impact is enhanced when computing the weak skin. For 
example, in the case of 48 Ca the spin-orbit enhancement amounts to about 0.03 fm; this is 
commensurate with the projected accuracy of the C-REX experiment. The impact of the 
spin-orbit contribution decreases with increasing mass number, so its effect in 208 Pb is small 
(0.006 fm) and well below the proposed ±0.05 fm accuracy. The spin-orbit contribution is 
in general small due to the strong cancellation between spin-orbit partners and between 
neutron-proton pairs in identical single-particle orbits. Given that it lacks the coherence 
displayed by the dominant (vector) contribution, spin-orbit currents have the largest impact 
on light neutron-rich (or proton-rich) nuclei with unpaired spin-orbit partners, as in the case 
of the H7/2 neutron orbital in 48 Ca. We note that in the RMF case where the cancellation 
among spin-orbit partners is incomplete, there is a further enhancement (of about 30%) 
relative to the non-relativistic predictions due to the presence of a strongly attractive scalar 
potential (the so-called "M*-effect"). 

In summary, accurately calibrated relativistic mean field models have been used to com- 
pute the impact of spin-orbit currents on the electroweak skin of a variety of nuclei. Given 
that spin-orbit contributions to both the charge and weak-charge radii may be significant 
and often as large as existent or projected experimental uncertainties, spin-orbit currents 
should be routinely incorporated into future calculations of both charge and weak-charge 
form factors. 
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Appendix: Single-nucleon electroweak currents 



In terms of the underlying quark vector currents, the electromagnetic and weak neutral 
currents displayed in Eq. Q are given by the following expressions [25] : 

4 = Yl Q^Qi = \u^u - l -d^d - U-fs , (A.la) 
f=i 

3 

■So = E atia^i = g>Yu + &*fd + glsYs , (A.lb) 
f=i 

where g w is the weak-vector charge of quark "f" expressed in terms of its weak isospin and 
the weak mixing angle. That is, 



A = 2Ti - 4Q f sin 2 9 W = { +1 " j S ^ *" * + °- 384 if f= | U ' C ^ (A.2) 

> -1 + | sin 2 6 W ~ -0.692 iff={d,s,b}. 



Note that we are assuming that the heavy-quark (c,b,t) content of the nucleon is negligible. 
Also note that the weak-vector charge of the quarks are to a good approximation equal to 
the negative of the electromagnetic charge of its weak isospin partner. This is the main 
reason why the neutral Z°-boson is an excellent probe of the neutron density. Assuming 
isospin invariance, namely, that the up(down)-quark distribution in the proton equals the 
down(up)-quark distribution in the neutron and that the strange-quark distribution is equal 
in both, the proton and neutron electromagnetic currents may be written in the following 
way: 

■JUP) = (A\p) = fv? - \V£ - \vit , (A.3a) 
= H J&Jn) = \v$ - \v* - \v* , (A.3b) 

where Vjf, V£, and Vjf are matrix elements of the respective quark vector currents in 
the proton. From the above equations one can determine V£ and V£ in terms of the two 
electromagnetic currents and Vjf. That is, 

V£ = 2J£ M ( P ) + J£ u (n) + V» , (A.4a) 
VS = Jem(p) + 2J^(n) + V* . (A.4b) 

In turn, these relations may be used to express the matrix elements of the weak neutral 
current in terms of the corresponding matrix elements of the electromagnetic current plus 
the strange-quark contribution. Inserting these relations into Eq. (A.lb) one obtains: 

J£ C (P) = 9UUp) + 3v n 4W + , (A.5a) 

J£c(n) = 9 n AiP) + 9lJUn) + & 0) V? , (A.5b) 

where proton, neutron, and singlet weak-vector charges are given — including radiative 
corrections — by [32J , 

gl = 2< + ^ = (1 - 4 sin 2 W )(1 + R p v ) « 0.0712 , (A.6a) 

9^ = <7v + H = -(1 + K) « -0-9877 , (A.6b) 

£? } = £ + 0v + 9l = -(1 + 4 0) ) « -0-9877 . (A.6c) 
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Note that these values are very close to the ones used in Ref. [2] which are the ones that have 
been adopted here (i.e., g% = 0.0721 and (7^ = —0.9878). Finally, electric and magnetic Sachs 
form factors for the weak neutral current may be expressed in terms of the corresponding 
electromagnetic form factors plus the strange-quark contribution as follows: 

GiuiQ 2 ) = 9 P v Gi M (Q 2 ) + 9*G£ M (Q 2 ) + 6 0) Gl M (Q 2 ) , (A.7a) 
G£ M (Q 2 ) = 9*G* M (Q 2 ) + 9 p v Gl M (Q 2 ) + ^GiuiQ 2 ) . (A.Tb) 

In particular, using these relations we obtain the following expressions for the nucleon 
weak-charge radii: 

9^1 = 9 P v r 2 p + 9>l + 6 0) r 2 s , (A.8a) 
9^1 = 9>l + 9 P v r 2 n + 6° ] r 2 • (A.8b) 

Similar expressions follow in the case of the nucleon weak magnetic moments: 

fi P = 9 p Vp + 9yVn + £,i 0) Vs , (A.9a) 
K = 9> P + 9>n + £i 0) l*s ■ (A.9b) 

Note that the mean-square strange radius and strange magnetic moment are defined as 

r 2 = -6 



dQ* 



and n s = GUQ* = 0). (A. 10) 

Q^=0 
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